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In this paper we provide strong evidence that there is no ambiguity in the choice of the horizon 
function underlying the Gribov-Zwanziger action. We show that there is only one correct possibility 
which is determined by the requirement of multiplicative renormalizability. As a consequence, this 
means that relations derived from other horizon functions cannot be given a consistent interpretation 
in terms of a local and renormalizable quantum field theory. In addition, we also discuss that the 
Kugo-Ojima functions u{p^) and w{p^) can only be defined after renormalization of the underlying 
Green function(s). 
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I. INTRODUCTION 

In 1977, Gribov [1] showed, in a saddle point approx- 
imation, that the restriction of the Euclidean functional 
integral to the Gribov region fl has far reaching implica- 
tions for the infrared behavior of the ghost and the gluon 
propagator. We recall that the region fl is defined as 
the set of field configurations fulfilling the Landau gauge 
condition and for which the Faddeev-Popov operator, 

M""' = -d^id^S-'+gr^'Al) , (1) 

is strictly positive. Therefore, 

n = {A"^, df,Al = 0, M"'' > 0} . (2) 

Later on, in a series of works, Zwanziger [2-5] elaborated 
on Gribov's approximation, being able to extend the 
previous results order by order at the quantum level. 
This resulted in an improvement of the Faddeev-Popov 
action which is now called the Gribov-Zwanziger action. 
In particular, the Gribov-Zwanziger action leads to a 
ghost propagator which is enhanced in the infrared 
region, a feature which has been confirmed by explicit 
two loop calculations in [6, 7], which constitute a non- 
trivial check of the predictions of the Gribov-Zwanziger 
formalism. 

Recently, it has been claimed [8-10] that the Gribov- 
Zwanziger action is plagued by a certain ambiguity. 
Depending on the choice of the so called horizon function 
[2], different results for the ghost propagator might be 
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found, namely: an enhanced or a non-enhanced ghost, 
according to [8-10]. Also, the recent lattice results seem 
to point towards a non-enhanced ghost[42] [11-15]. A 
natural conclusion would seem to be that one should take 
the horizon function which leads to the non-enhanced 
ghost. However, this does not agree with the original 
results by Gribov and Zwanziger and therefore asks for 
an explanation. 

In this paper, we would like to clear the situation. 
We shall show that there is no ambiguity in the choice 
of the horizon function. The correct form of the horizon 
function is the one originally constructed by Zwanziger 
[4], and is clearly dictated by the renormalization 
properties of the Gribov-Zwanziger action. We stress 
that renormalization is of paramount importance for 
defining meaningful Green functions. A dynamical 
improvement of this Gribov-Zwanziger action, consistent 
with the renormalization, consequently allows to obtain 
the non-enhanced ghost, as discussed in previous work 
[21, 22], giving results compatible with other analytical 
approaches [18, 23-26], based on the Schwinger-Dyson 
formalism, see also [27]. 

As a corollary of the present analysis, we shall elaborate 
on the meaning of the Kugo-Ojima functions u(p^) and 
defined as follows [28, 29] 

J d^xd'ye^'^^^-y^ {igfatcA'^c'){x){gfaMA';c'){y))^^^ 
= g,Mp') + ^w{p^) (3) 

and their meaning at the level of renormalization, shed- 
ding more light on certain claims in [8-10]. We shall also 
discuss about several results obtained in the literature 
[8-10], where the other choice of the horizon function 
was investigated. As we shall show that this particular 
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choice cannot be consistently introduced within a local 
and renormalizable framework, it follows that its con- 
nection with the ghost propagator's behaviour, encoded 
in the quantities u{p'^) and w^p^), cannot be properly 
defined at the quantum level. 

The paper is organized as follows. In section II, 
after giving a detailed account of the aforementioned 
claims [8-10] about the existence of a possible ambi- 
guity in the choice of the horizon term, we shall show 
that there is only one possible horizon function which 
is dictated by the renormalization properties of the 
Gribov-Zwanziger action. This shall be extensively dis- 
cussed in sections III and IV, thereby showing that there 
is in fact no such ambiguity. In section V, we shall spend 
some words on the Kugo-Ojima Green function (3) and 
on an alternative one which, in contrast with expression 
(3), turns out to be renormalizable and closely related 
to the horizon function. We shall end this paper with 
a discussion, by emphasizing our possible explanation 
for the recent lattice results. In the Appendix A, we 
have provided the complete proof of the renormalization 
of the Gribov-Zwanziger action, needed for the analysis 
performed in sections III-V, whereby we have also taken 
the opportunity to correct a small mistake concerning 
a minor statement in our previous works [22, 30, 31]. 
However, none of our results [21, 22, 31] are affected by 
this. 



II. SURVEY OF THE ISSUE CONCERNING A 
POSSIBLE AMBIGUITY IN THE CHOICE OF 
THE HORIZON FUNCTION 

In [4] , it has been shown that the restriction to the Gri- 
bov region f2 can be established by adding the following 
nonlocal term 

5h = y d'^xh{x), (4) 

to the standard Yang-Mills action. 

Sym + -Sgf , (5) 

with S'ym the classical Yang-Mills action and S^i the Lan- 
dau gauge fixing 

Sym = ^'''^^Ili^P^i^ ' 

Sgf = J d'x {b'^d^A;+rd^D;''c') . (6) 

The nonlocal expression h{x) is called the horizon func- 
tion. In [8] , it has been argued that there are in fact two 
possible choices for the horizon functions, namely 

h,{x) =7^1 dVr'''A'^{x){M-'r\x,y)r^Al{y) , 

(7) 



or 

h2{x) = d''yD';,''{x){M-'r''{x,y)D'';{y), (8) 

with A4 the Faddeev-Popov operator 

Al""^ = -d^D'^f, (9) 

and 

Df = d^S'^^' + gr^^Al (10) 

the covariant derivative. The parameter 7 is not free, 
but is fixed by a gap equation, known as the horizon 
condition, which reads 

{h{x)) = d{N^ - 1) , (11) 

with d the number of space-time dimensions. This condi- 
tion ensures that the domain of integration in the Feyn- 
man path integral has been restricted to field configura- 
tions belonging to the Gribov horizon fl. In this paper, 
we shall motivate that there is only one possible hori- 
zon function, the one which was originally intended by 
Zwanziger[43] [4], given by 

5h = lini J d'^xh^iix) = \iinj d'^x J d'^y 
x{D';,%xh^ix))iM-'r\x,y){D'';{y)j'iy)) , (12) 
whereby 7(2;) is defined through 

^^{z) = e'^^72 . (13) 

The limg^o operation corresponds to replacing the 

space time dependent 7^(2:) with the constant Gribov 
parameter 7^. We observe that this horizon function 
shares a great resemblance with hi{x) and h2{x). It 
is worth to point out here that the limit, lime^o, in 
expression (12) is meant to be taken after an appropriate 
localization of the horizon function, a point which we 
shall outline in detail in what follows. 

Let us try to explain here how it is possible that 
some doubts have arisen about which horizon function 
is the correct one. Firstly, one can easily see that 
setting 7(x) in /i3(x) immediately equal to a constant 
parameter, "f{x) = 7, which agrees with switching the 
limit limg^o a-nd the integration signs, we obtain the 
horizon function hi{x). One can therefore appreciate 
that the difference between hi{x) and h3{x) is very 
subtle. Comparing h2{x) with h-T,{x), it is apparent that 
the space time dependent function 7^ (y) has been pulled 
out of the horizon function, an operation which is not 
completely justified, due to the nonlocal character of 
the kernel {M~^)°'^{x,y). A second source of a possible 
ambiguity can also arise when we try to localize the 
horizon function hz{x). Indeed, we can localize the 
action S'ym + Sgf -|- Sh by introducing some extra fields. 
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Looking at the following standard formula for Gaussian 
integration for bosonic fields 

Cdet^-^exp J dVy J^{x){A-'r\x,y)4{y) 

+ 1 d^'x {^'^j;{x) + ip%x)J^{x))] , (14) 

we observe that we can get rid of the inverse of the 
Faddeev- Popov operator in /i3(a;) by introducing new 



fields. For every index i, defined by . . = . . we 



can write for h3{x) 

expi- d'^xh3{x)\= Yl det{-M) [d(p][d<f] 

exp|^lim[y d'^x J d''yTptix)M"\x,y)cp'i{y) 

+ J d'^x (D^ixh^x)) ^^{x) + (D^ixh^x)) Jpfix) 

(15) 

whereby we have introduced a pair of complex conju- 
gate bosonic fields {'<^'^,f1f) — {'<ff,(pf). We can then 
also lift the determinants det{—A4) into the exponential 
by introducing a pair of Grassmann fields (uJ'^'^ , u'^'^) = 
(ZJfjWf). Making use of the standard Gaussian formula 
for Grassmann variables 

C(det A)exp J d'^xd'^y J^(x)(^-i)»''(x, 2/)4(y) 

= y"[dw][daJ]exp j d'^xd'^yuf{x)A''''uj\y) 



+ / d^x (J»(x)a;-(x)+Z:;"(x)J^(x)) 



(16) 



whereby we set the sources and equal to zero, we 
obtain 

exp i^- J d''xh3{x)j = n y [^^] [^'^j [^^] [^^] 

exp[y d-^x J d''y{lpt{x)M"\x,y)ip'dy) 

^utix)M^Hx,y)co',{y)) 
+ \\mj d^x {Dt{x)^\x)) ^1{x) 

+ (D«(x)7^(x))^«(x)" . 
The new localized action thus becomes 

Sqz = S'q + S'^ , (17) 



with 

S'o = 5ym + % + J d'^x {r;d.Dfvl' - uj^'^d.Dfuj';) , 

(18) 

and with 

^7 = -1™/ d'x[{D;^{x)^\x))^l''{x) 

+ (D-(x)7^(x))^-(x)] 
= limyd'^X7'(x)i?-(^-(x)+^-(x)) 

= 7'/d'^a;I?-«(x)+^-(x)). (19) 

Notice that, as already remarked, the limit 6 ^ 0, in 
equation (19) has been performed after localization. As 
one can see from (13). taking this limit is equivalent with 
setting 7^(x) equal to the constant 7^. As at the level 
of the action, total derivatives are always neglected, S!y 
becomes 
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Jd'^xgr^-AK^l^+lpl-) . (20) 



From this expression we can easily understand why 
certain doubts have arisen. Starting from the first 
horizon function hi{x) given in (7) and undertaking the 
same procedure, we would end up with exactly the same 
action S!y. This can be understood as we have neglected 
the total derivatives. In conclusion, although the local 
actions derived from hi{x) and h^ix) are the same, at 
the nonlocal level they are clearly different, which also 
follows from [8]. This is important when one is doing 
manipulations at the level of the nonlocal action as has 
been done in [8]. 

Let us now translate the nonlocal horizon condi- 
tion (11) into a local version [4]. The local action 5gz 
and the nonlocal action Sym + Sg{ + Sii are related as 
follows, 



/ 



[d^] [d6] [dc] [dc]e-('5^"+'^^f +'^'-) 

= / [dA][d6][dc][dc][d<y?][d^][da;][diU]e-'^°^ 



Next, we take the partial derivative of both sides with 
respect to 7^ so we obtain, 

-2j'{h) = {gr'^A^i^"; , (21) 

for both horizon functions hi{x) and /i3(x). We recall 
that (9^99"") = and {dfjp"'"') — 0, meaning that both 
horizon functions hi and give rise to the same local 
horizon condition. Using these expressions and assuming 
that 7 ^ 0, we can rewrite the horizon condition (11) 

{gr''A';^{^l^ + 1;)) + 2^H{N^ - 1) = . (22) 
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By adding the vacuum term 

to S!y, we can write the horizon condition as 

, 



dT 



(24) 



with r the quantum action defined as 

e"^ = j [d^je"^'^^ , (25) 

where /[d$] stands for the integration over all the fields. 

For the Gribov-Zwanziger action to be renormaliz- 
able, it necessary to perform a shift over the field wf, 

see [4], 

+ J d''z{M-'r\x,z)gUked^[D';^c%z)ip\z)] , (26) 

so that the action becomes 

Sgt. = So + , (27) 
whereby Sq has been replaced by 5*0 

So = S'o + J d^x , (28) 



and the vacuum term is now included in S-y 



S^ = S'^ + J d'^x jUiN^ - 1) 



(29) 



The renormalizability of ^qz has been proven in the Ap- 
pendix A. We would like to stress that this is far from 

being trivial, especially since no new parameter is needed 
to take into account the divergences of the vacuum term 
we introduced by hand in equation (23). In addition, the 
algebraic formalism employed in the Appendix A also 
gives a more clean argument why the extra term appear- 
ing in equation (28) is necessary, without the need of 
performing the nonlocal shift (26). 



III. THE COMPOSITE OPERATORS, PART I 

We shall provide a strong argument which illustrates 
that only the horizon function /is (a;) possesses a clear 
meaning at the quantum level. For this, we first need to 
demonstrate the following equality 

{f{A){M-'r\x,y)) = {f{A)c'^{x)c\y)) 

= -^{f{A)ujf{xp'M)m 



where f{A) stands for an arbitrary quantity depending 
on the gauge fields A^, and is the number of colors. 
Let us start with 

(/(^)c"(x)c*(y)) = y'[d$]/(A)c«(x)c''(j/)e-^°M-31) 

We can rewrite this expression by adding corresponding 
sources for the ghost and the gluon fields, 

{f{A)c'^{x)c'iy)) 



f 

with 

'S'so 

= / d-'a; 



\SJaJ 84{y) 5J^{x) 



[d$]e 



-'S'GZ + 'S'f 



sources 



sources=0 



(32) 



J- 



J. 



(33) 

and with the full expression for 5gz given in (27). We 
can perform the integration over the ghosts c, c, w, w as 
this is just a Gaussian integration. The relevant piece of 
the action is given by 

At"" gfaudMDf) 
M"^ 



Sgz = J d'^x [ujf c«] 









c''_ 



Kab 



Making use of the formula (16) gives 
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(34) 



(wa) 54{y) 5J-{x) 

\—b 

[d$]exp<j [J-^ 4] {K-'r" 



{f{A)c'^{xrc\y))f 

j 



sources=0 



(35) 



The matrix {K ^)"° can be computed as 



-l\ab 







X 



(36) 



whereby x is some function of the relevant fields. From 
(32), we now see that 

</(/l)c°(rt)c'(!,)> = jmS(A){M-^)-^e-'"- (37) 
(;,^K„^,„,) = _^/(J_) * _L_ 



while 



y[d$]exp|[j-^ J-]{K-^r' 





-I 







sources=0 



= -:^{f[A)ujt{x)^\{y)) 
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which proves relation (30). 



A. The horizon function hs 

We shall now prove that the expectation value (/is (a;)), 
see equation (12), is renormalizable. We can apply the 
formula (30) to the horizon condition {hs^x)), yielding 

j d'^x {h3{x)) = lim j d'^x j d'^y 

{{Dl%x)^\x)) {M-'r\x,y) {Dl%y)^\y))) 

= -\imJ d'^x J d'^y 

{{D^lfixh^ix)) utixpliy) {Dl^iyh'iy))) 

= -Jf^' j / <^''y{{D,u:inx){D,Winy)) . 

(38) 

In order to check whether this horizon term {x) is well 
defined, the correlator {{D^uji)°'{x){D^uJi)°-{y)) includ- 
ing the composite operators {Df^uJi)°'{x) and {Df^uji)'^{y) 
should be renormalizable. In fact, this turns out to be 
the shown in the Appendix A, where a detailed 

account of the algebraic renormalization of the Gribov- 
Zwanziger action has been provided. Moreover, we have 
also proven that the action Eqz (see equation (A5)) is 
renormalizable. From this, one can immediately obtain 
the correlator {{D^ujiY{x){D^'UjiY{y)) by deriving the 
action Eqz (A5) with respect to the sources[44] N^{x) 
and U^iy), 



[d$] 



S 



SN-ix) 5U^{y) 



all sources = 



= j m [-9r''\D^c)\x)^i{x) + {Dujinx)] 



{DUir{y)e-''^- 
= {-gr'%D^c)\x)^U^)iDZJ,r{y)) 

+ {{Duj.nx)iDZJiny)) . (39) 

We shall now show that the first correlator in expression 

(39) vanishes, 



{{9r'''{D,cf^m^){Duir{y)) = 



(40) 



In fact, equation (40) belongs to a more general class of 
Green functions which are all zero, namely 



(e(x)A(2/)) = 



(41) 



with Q{x) a function of fields not containing the field u, 
while A{y) is a function containing oJ. We can show that 

all these Green functions are zero by using an elementary 
diagrammatical argument. It is impossible to construct 
any diagram which has an oJ leg starting from a space 
time point y which has to be connected in some way to 
a space time point x, where no w leg is present. Indeed, 



every ZJ requires an uj leg to propagate, w in his turn 
shall always produce another lu leg in all vertices as 
can be seen from the action (27). Moreover, the field uJ 
needs again another u) leg to propagate. Therefore, an 
u leg is required in the space time point y to close the 
diagram. As a consequence, all Green functions of the 
type of equation (41) are zero. 

Hence, we conclude that the Green function 
{{DojiY{x){DWiY{y)) is multiplicatively renormalizable 
as follows from 



{{Du:ir{x){Duj,Y{y))o = Zu'Z-N {{DiOiY{x){DUiY{y)) , 
whereby 



y—ly—l y'-Z-^y 7—1 



_ -71/2 

r ■^jv 

see expression (A41). 



(42) 
(43) 



B. The horizon function h\ 

We shall now prove that the first horizon function (7) 
implies a horizon condition which is not multiplicatively 
renormalizable. In an analogous fashion as in the previ- 
ous subsection, we can write 



/ 



d^xl^hx{x)) = 



' j^'xj dS{9T'"'Al{x){M-^Y\x,y)f"'Kiy)) 

~j'J d'x J d%{{gr''^A';u^!^){x){gf''AloJt){y)) . 



(44) 



We can demonstrate that the composite operators 
gfo-kcj^k^a g^j^i^ gf^^'^A^^uj\ are not renormalizable. This 
is due to the fact that, at the quantum level, those com- 
posite operators will unavoidably mix with the operators 



d,.. 



which have the same quantum numbers. 



For this, we need to consider the operators c^^wf and 
gfakbA^^<jj\ as separate operators, each coupled to their 
own source, instead of D'^ijj\ being coupled to the single 
source U^^, and similarly for 9^aJ" and g fakbA'^oJ^ . In 
Appendix B, we have given an alternative proof of the 
renormalizability of the Gribov-Zwanziger action, where 
we have coupled the following sources to the following 
operators, see equation (B2), 



gfakbj^k^b_gf^^^j^ba^a^c 



U" 

u',: 



gr^'A^^cot 



(45) 



In (53), we have foimd that the sources U and U' mix, 
as well as N mixes with N'. Taking the inverse of this 
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matrix yields 



" u ' 




ail 




'Uo' 


U' 










' N ' 




' Z-^ ai 




'No' 


N' 




z'J\ 







(46) 



Wc recall that insertions of an operator can be obtained 
by taking derivatives of the generating functional 
Z''{U, U', N, N') w.r.t. to the appropriate source. For 
example, 



{gf A^u,)o 

5N 5Z%U, U', N, N') SN' 6Z%U, U', N, N') 



+ 



5N 5N^ SN 

so that 

(5r'=''4cU?)o = ai(a^iI;«) + 4^'(ff/"''4^i)-(47) 
We can do the same for the other operator 

= a.id^col) + Zgigr'^'Ay, ^ gU,D^;^<^>^) (48) 

Taking the partial derivatives of given in equation 
(Bl) w.r.t. iV;"(x) and U'^(y), we find 



/ 



m 



all sources — 



5N'^{x) 5U'^{y) 
= J m {gr'-'Ay, - gfabcDl^c^^t) (^) 

Agr''Aluj1){y)e-^^- 
= {{gr''Ay,){x){gr'^AlzJt){y)) 

- ((5/„6e<cVD(^)(5/'^'^4^-)(2/)> • 

The correlator {{gfabcD'lfc'^ip^){x){gP'''>A''^uj''^){y)) also 
belongs to the class (41) and is therefore equal to zero. 
However, the remaining Green function 
(^{^g fo-'^b A^^u)\){x){g f"-^'^ A^^uj1){y)) is no longer multiplica- 
tively renormalizable. Indeed, due to the mixing we have 
found, the bare correlator can be written as follows 

{{gr^'A';u:t){x){gr"^Alu;'i){y))^ 

= a.z'f {{d,uj1){x){gr'^A'^uj1){y)) 
+a^Zg{{gr'''Ay,)ix){d,uJt)iy)) 

+Z,Z\/'{{gr'''>Alojt){x){gr'^Aluf,){y)) . (49) 

This is a strong argument why one should not rely on 
this first horizon hi function. When renormalizing this 

Green function, the "missing" terms stemming from the 
covariant derivative re-enter again. 



IV. THE COMPOSITE OPERATORS, PART II 

A. The horizon function ha 

Let us now look at the composite operators which are 
involved in the horizon term h^ after having performed 

the necessary localization. In particular, wc are inter- 
ested in the composite operators appearing in expres- 
sion (19) before setting 7(x) equal to a constant, namely 
D'^^ip'^^ and D^^Tp''^. If the horizon function is well de- 
fined, we expect these composite operators to be renor- 
malizable. We can check this again from Appendix A. 
Firstly, from (A5) we can deduce 

(z)«Vf)o = z]^\D;'^t) = zy'z]/\Df^t).m 

Secondly, we can do something analogous for Dflfl We 
see that this operator is a linear combination of two 
renormalizable composite operators, namely the compos- 
ite operator coupled to i.e. -Dfjp'l + gfabcD^^d^tdl 
and to R^, i.e. —gfabcD'^c'^cJi. Luckily, the two compos- 
ite operators have the same renormalization constant, as 
without this property, the linear combination would not 
be renormalizable. Therefore, 

= Zy\DfTp\), (51) 

A/2 7I/4 



whereby Zy^ = Zg^Z^J^'^ as can be foimd in (A41). In 
conclusion, bo 
are renormalizable. 



conclusion, both composite operators D'^^'ip''* and D'^^'ip''^ 



B. The horizon function hi 

When taking the limit 9{x) 0, we can drop the total 
derivatives. However, the remaining composite operators 
gfakbA'^ip'l and gfakbA'^Tp'^^ are not multiplicative renor- 
malizable. We can prove this in an analogous fashion as 
in the previous section. In Appendix B, we have con- 
sidered the operators 9^<y3f and gfakbA'^ip^ as separate 
operators instead of being coupled only to one source, 
and similarly for d/jTpf and gfakbA^Tp\. To that purpose, 
in the Appendix B, we have coupled the following sources 
to the following operators, see equation (B2): 

{gfakbAlTp", - gUcDl^Cufi) ^ F;-, (52) 

and we have reworked out the complete algebraic renor- 
malization in this alternative setting. In (53), wc have 
found that the sources M and M' mix, and also V mixes 
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with v. Taking the inverse of this matrix yields 



" M ■ 








'Mo' 


M' 




^r^y\ 






' V ' 




" ^-l/2„-l/4 




>o" 


V 




z^zT 




To. 



(53) 



As in the previous section we have that 



{gr''A';^% = a,{d,^t) + Z-'I^Z~'/\gr'^^Al^\) . 

(54) 

We can do the same for the other operator. Only some 
care has to be taken as the source V^"' couples to a sum of 
two operators [gfakbApp\-gfabcD^ifc'^'^'i)- Wc therefore 
have to subtract the operator coupled to the source R'^, 
namely 

{gfakbAli^)o 

= {gfakbAy, - gU^Dl^c^ujDo - {-gfabcOl^c^ufl)^ 
= a,id,^) + Z-^^l^Z'^^l\gr^''Alip\ - y/abc^C^SJ- ) 

-zy^z)i\~g!ab.D\'c<'^;) 

= a,(d,^) + Z;'/'Z-'/\gr'''>Al^1) . (55) 

We can thus conclude that the operators gfakbA''ip\ and 
gfakbA^Tp\ are not multiplicatively renormalizable and 
mix with the operators dfj_Lp^ and dfjjpf, respectively. 
Therefore, one should keep in mind that the limit 9—^0 
has to be taken as the final step, and that it can only be 
taken at the local level. Furthermore, the mixing we have 
found tells us that one should always leave the covariant 
derivative of a field "in one piece". As a consequence, 
much care has to be taken at the nonlocal level when de- 
riving all kinds of results, as we shall now explain in the 
next section. 



V. THE KUGO-OJIMA GREEN FUNCTION(S) 
AND THE LINK WITH THE HORIZON 
CONDITION 

As explained in [33], there exists a close link between 
the Gribov-Zwanziger formalism and the Kugo-Ojima 
analysis of gauge theories. An important role in the 
Kugo-Ojima work is played by the parameter «(0), which 
is the zero momentum limit of the function u{p^), which 
can be extracted from 



d'^xe'f^ {D''^'^c'^{x)Dl''r{0)) 



be—ef 



= 5' 



ab 



Shu - 



Pt^Pu 



(56) 



In order to have a finite quantity u{p'^), the previous 

Green function should of course be renormalizable. In 
the case of pure Yang- Mills theories, one can indeed prove 



this, but this was achieved by also using anti-BRST in- 
variance. We refer to the literature [34, 35] for details. 
In the presence of the restriction to the Gribov region 
il, we even lack the concept of the anti-BRST symmetry 
as far as we know. We shall therefore present a slightly 
different argument here. Using the correspondence (30), 
we may also write 



<p1 



PixPv 
p2 



(57) 



from which the renormalizability immediately follows, see 
equation (42). If one wishes to introduce bare quantities, 
one can write down 

j A^xe''^^ ipfu^i\x)Dl^uf^m^ 

u{p') 



= 6"'' 

Z-' 



PtiPv 



P^lPu 



= z-H"^ 



d'^xe'P" {Dfoji\x)Dl''uJ''^{Q)) 



Pf^Pu 
p2 



p2 



,(58) 



but this is just a formal relation. It nevertheless 
teaches us that one cannot simply speak about the 
(multiplicative) renormalization of u{p^) in terms of a 
bare Kugo-Ojima function uo(p^). 

As a corollary, using the identification (38), we de- 
rive the following correspondence between the horizon 
function and the Kugo-Ojima parameter u(0). 



-{h,) = iN^-l)[{d~l)u{0)-l] 



(59) 



which hereby rigourously proves the relation which was 
one of the starting points of the paper [33] . 

In [29], it was shown that one can parametrize the 
ghost propagator, defined as follows 



{c^{-p)c\p)) = 6-^G{p^) , 
in terms of[45] 



p^(l + u{p^) + w{p'^)) 



(60) 



(61) 



This relation was also discussed in [4, 8-10, 26, 36]. Us- 
ing this relation, one can again derive a formal correspon- 
dence between bare and finite quantities, 

[1 + w(p2) + w{p'')]a = Z-\\ + w(p2) + w{p^)\ . (62) 

but as already noticed in [36], it does not allow to 
separately speak about (1 -|- u(p^))o and w^i^^. 
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In [8], it was claimed that the funetion w(p^) en- 
joyed the property of being not-renormalized. We are 
unable to understand this claim, especially since we do 
not know what the corresponding bare quantity wo{p^) 
would be. Let us explain. Due to expression (30), we 
can rewrite (3) as 



- 1 ^d'^ye'f (-^) {{gU,Alu'^){x){gUeAtTni){y)) 



ipi 



(63) 



We see that this correlator is, up to the Lorentz structure, 
equal to the one investigated in equation (49). Also this 
correlator is not multiplicatively renormalizable, as we 
find 

{{gr'''Aluj'^){x){gr'^AiuJmy)), 
= a.Z'J' {{d,u:'}){x){gr'^Aiz3-i){y)) 
+a,Z,{{gr'''Alu;'>)ix){d,i::t)iy)) 
+ZXJ' {{gr'''Ay,){x){gr'^Aiuf,){y)) .(64) 

It is therefore of no use to define m(p^) and w{p'^) at the 
bare level, as we cannot introdiicc their corresponding 
bare counterparts. We can conclude that the definition 
of w(p^) and w(p^) should be considered as a definition 
at the level of renormalized Green functions, in which 
case the parametrizations (56) and (61) make perfect 
sense. Due to the mixing, we do not see how we can 
maintain the definition (63) at the renormalized level. 
In our opinion, one must first define u{p'^) through 
equation (56), which then allows to extract a value for 
w(p^) from the ghost propagator via the relation (61). 

We end by commenting on the results of [8 -10] 
which gave a connection between the horizon function 
hi given in equation (7) and the functions u{p^) and 
w{p'^). Using the nonlocal expression (7), one manages 
to write down a quite complicated connection. However, 
as we have clearly shown in this work, it is impossible 
to treat hi at the quantum level, as quantum effects 
make necessary the presence of the derivative operators, 
see the previous section. At the end, only the horizon 
function h^ given in equation (12) is meaningful, and as 
such only the connection (59) survives. In contrast with 
the outcome of [46] [8], which turns out to be ill-defined, 
imposing the horizon condition (11) would lead to the 
infrared enhanced ghost as follows from the results (59), 
(61), and using the fact that w{0) = 0. This has been 
confirmed explicitly in [6] to two loop order, and it has 
also been found using a Schwinger-Dyson approach in 
[26]. 

The question that remains then is whether the 
ghost propagator is necessarily infrared enhanced when 

enforcing the horizon condition? Indeed, starting from 
the standard Gribov-Zwanziger action, one shall find 



an enhanced ghost propagator. As already indicated in 
the introduction, this does not seem to be in agreement 
with recent lattice results. However, after one takes into 
account the nontrivial dynamics of the Gribov-Zwanziger 
action by including the possibility of having nonvan- 
ishing condensates, one is able to find different results 
[21, 22, 38]. By including for example the dimension 
2 condensate {iffipf —uJfujf), a non-enhanced ghost 
propagator is found within the so called refined Gribov- 
Zwanziger (RGZ) formalism [21, 22, 38]. Simultaneously, 
the same condensate also leads to a gluon propagator 
that does not vanish at zero momentum. Moreover, we 
mention that the form of the gluon propagator stemming 
from the refined Gribov-Zwanziger framework enables us 
to fit the lattice data in a very reasonable way [38]. We 
notice that the works [8-10] are unable to say anything 
about the gluon sector. 



VI. CONCLUSION 

The main conclusion of this paper is that the renor- 
malizability requirement of the Green functions of a 
quantum field theory is of paramount importance. It 
dictates whether the Green function one is considering 
makes sense or not. Once having obtained a finite, 
renormalized Green function, one can address the 
question of how the renormalized Green function can be 
parametrized in terms of certain quantities. It is of no 
use, or even incorrect, to discuss the renormalization of 
such parameters directly. 

An important ingredient to discuss and interpret 
the underlying renormalizable quantum field theory is 
the issue of locality. If one wishes to study nonlocal 
theories, a prerequisite appears to be that one can write 
down a local version of the theory, to which the usual 
tools of local quantum field theory apply [39]. Once 
having established the necessary features of a quantum 
field theory at the local level, one can try to go back and 
see what the corresponding nonlocal version looks like. 

In particular in this paper we have first shown that 
the Gribov-Zwanziger action is unambiguous. There is 
only one possible horizon function at the nonlocal level, 
which is given by 

J d'^xhix) = Vim J d'^xhaix) = \imj 

(D-(.x)7^(x)) {M-'rix,y)iD^^{yh\y)) . (65) 

with 7^(2) = e*^^7^. We have motivated in two different 
ways that this is the only correct horizon function to 
start any discussion from. Our analysis was based on 
the renormalization properties of the corresponding local 
version of the operators involved. 

In addition, we have also elaborated on the defini- 
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tions of the Kugo-Ojima parameters u{p'^) and w{p'^), 
and we have argued that it only makes sense to talk 
about these parameters at the level of renormalized 
Green functions. We reanalyzed the connection between 
u{p^) and the horizon function, and motivated that the 
relation (59) is the only correct one, in contrast with 
other results made available recently [8]. We believe 
that any relation derived from other horizon functions 
should be looked upon with great caution [8]. 
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Appendix A: Renormalization of the 
Gribov-Zwanziger action, option 1 

In this section we shall repeat the proof of the renor- 
malizability of the Gribov-Zwanziger action. We shall 
also have the opportunity to correct a minor mistake 
present in previous works [22, 30, 31]. This is the rea- 
son why we shall repeat the full algebraic renormaliza- 
tion in all details. In particular, our current analysis 
also presents a slightly improved version of the original 
renormalization proof [4, 40]. However, let us stress that 
none of our results are affected by this minor error, like 
e.g. those presented in [22, 31]. 

1. The stcirting action and the BRST 

We start with the following action, 

Sgz = So + S^, (Al) 

with 

So = Sym + 5gf + J d^x {iptd^ {Df^\) 

+r^^^»^^^ + -(7V2-l) 7')(A2) 



First, notice that we have simplified the notation of the 
additional fields (^^'^j i^^^, w^*^, w^'^) in Sq as 5*0 displays a 
symmetry with respect to the composite index i = {fi, c). 
Therefore, we have set 

{^fi 'V>n ^^ix ^^ti ) = VPi,Vi,i^i,^i) ■ (A3) 

The conventional BRST variations of all the fields are 
given by. 









= l9f 






sb" 


= 0, 




— ^i, ' 




= 0, 




= 'Pi , 




= 0. 



One can check that the BRST symmetry is softly broken 
for the Gribov-Zwanziger action [2, 22], 

sSgz = s{So + Sj) = s{S^) 

= 57' / d^xr'^ {Ay; - {Dl^c^) {if; + f;)) . 

We notice that the breaking is due to the 7 dependent 
term, S^. 

In order to discuss the renormalizability of 5gz, 
we should treat the breaking as a composite operator 
to be introduced into the action by means of a suitable 
set of external sources. This procedure can be; done in 
a BRST invariant way, by embedding S'gz into a larger 
action, namely 

^Gz = Sym + Sgi + So + Sg, (A5) 

whereby 

Ss = sj d'^x {-U^Dff, - V^'Dfufi - U;v;' 

= J d^x {-M^Dff, - gr^^U^Dl^d'ip'l 

I rrai Tkab, b Ajai r\ab—b Trai r\ab—b 

+gr'"'V/^' D'lfc:^ oJ^^ - M^V"' + U;n; 
+R;g.r''W';c''Zo^,+T;gUcDl'c''r^) ■ (A6) 

We have introduced 3 new doublets {U^, Mf), (V^"*, 
Nf) and {T;, R;) with the following BRST transfor- 
mations, and 

sU; = M;, sM; = 0, 

sT; = Rf , sR; = . (A7) 

We have therefore restored the broken BRST at the ex- 
pense of introducing new sources. However, we do not 
want to alter our original theory (18). Therefore, at the 
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TABLE I: Quantum numbers of the fields. 







c" 


c" 


6° 


^1 


^■ 






dimension 


1 





2 


2 


1 


1 


1 


1 


ghost number 





1 


-1 











1 


-1 


Q /-charge 














1 


-1 


1 


-1 



TABLE II: Quantum numbers of the sources. 



rjai 




prai 


yar 


r>ai 


rjiai 

-'m 






2 


2 


2 


2 


2 


2 


3 


4 


-1 





1 








-1 


-1 


-2 


-1 


-1 


1 


1 


1 


1 









end, we have to set the sources equal to the following 
values: 



I phys 



phys 



f^^ I phys 



yab\ 

I phys 



• I phys 



= 0. 



/^^ I phys 



It is exactly here that we have committed a mistake. In 
previous papers, we have forgotten to introduce the dou- 
blet (T", i?"), and therefore, the intended physical limit 
did not exactly reproduce the original action (Al). In the 
original article [4] , these two sources were also not intro- 
duced. When taking the physical limit, an extra term 
was generated, which was then removed by doing a (non- 
local) shift in the uj field. Here, wc have circumvented 
this unnecessary shift by introducing the doublet (T"^', 



1. The Slavnov- Taylor identity is given by 
<5(S'gz) = , 

with 



(A13) 



GZ 



GZ 



5^1 

+ R': 



GZ 



2. The U{f) invariancc is given by 



(A14) 



2. The Ward identities 

Following the procedure of algebraic renormalization 
[41], we should try to find all possible Ward identi- 
ties. Before doing this, in order to be able to write the 
Slavnov- Taylor identity, we first have to couple all non- 
linear BRST transformations to a new source. Looking 
at (A4), wc see that only A^^ and c° transform nonlin- 
ear ly under the BRST s. Therefore, we add the following 
term to the action Egz, 



with and L" two new sources which shall be put to 
zero at the end. 



VIphys ^ Iphys ~ " • 



(AlO) 



These sources are invariant under the BRST transforma- 
tion, 

sK^ = 0, sL" = . (All) 

The new action is therefore given by 

^GZ = ^GZ + Sext ■ (A12) 

The next step is now to find all the Ward identities 
obeyed by the action T,'q2,- We have enlisted all the iden- 
tities below: 



- M"^ 



^ 5R?,i ST'f 



By means of the diagonal operator Qf = Uu, the 
z-valued fields and sources turn out to possess an 
additional quantum number. One can find all quan- 
tum numbers in Table I and Table II. 



3. The Landau gauge condition reads 



6n 



GZ _ a Aa 



4. The antighost equation yields 



•^^GZ _|_ g _<^^GZ 



SK"; 



= 



(A15) 



(A16) 



5. The linearly broken local constraints yield 



AS' 



GZ 



O '^^GZ j_ „ P rpdi <^^GZ 



^^'gZ I Q '^^GZ 



gf 



abcr^b '^^GZ 



5¥ 



(A17) 
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6. The exact TZij symmetry reads 



(A18) 



with 



+V' 



" SRI' 



(A19) 



7. The integrated Ward identity is given by 



Here we should add that due to the presence of the 
sources T^* and R'^\ the ghost Ward identity [41] is bro- 
ken, and we are unable to restore this identity. For the 
standard Yang-Mills theory, this identity has the follow- 
ing form 



(A21) 



with 



and 



= /'^'K^^^^^'^''^) ' ^^^^^ 
Ki=9j d'^/"'^ (K^l - L^'c") , (A23) 



a linear breaking. However, it shall turn out that this is 
not a problem for the renormalization procedure being 
undertaken. 



with B the nilpotent linearized Slavnov- Taylor op- 
erator, 



B 



5K^ (5^« 



5c°- 5L° 



■ n: 



and 



B^ = 0. 



R 



S 



2. The U{f) invar iance reads 

3. The Landau gauge condition 



GZ 



0. 



4. The antighost equation 



'^^GZ , p) ^^GZ 

+ " SK- 



5. The linearly broken local constraints yield 



(A26) 



(A27) 



(A28) 



(A29) 



(A30) 



^GZ 



, (A31) 



3. The counterterm 

The next step in the algebraic renormalization is to 

translate all these symmetries into constraints on the 
counterterm Sq^i which is an integrated polynomial in 
the fields and sources of dimension four and with ghost 
number zero. The classical action Sq^ changes under 
quantum corrections according to 

^GZ ^ ^GZ 

+ /iS^z , (A24) 

whereby h is the perturbation parameter. Demanding 
that the perturbed action {T,Q^ + hT,Q.^) fulfills the same 
set of Ward identities obeyed by Sq2, see [41], it follows 
that the counterterm is constrained by the following 
identities. 

1. The linearized Slavnov- Taylor identity yields 

SS^z = , (A25) 



= . (A32) 

6. The exact TZij symmetry reads 

Te^E^z = , (A33) 
with TZij given in (A19). 

7. Finally, the integrated Ward identity becomes 

Now we can write down the most general counterterm 
of d = 4, which obeys the linearized Slavnov- Taylor 
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identity, has ghost number zero, and vanishing Qf num- 
ber, 



(A35) 

with ao, . . . ,aii arbitrary parameters. Now we can un- 
leash the constraints on the counterterm. Firstly, al- 
though the the ghost Ward identity (A21) is broken, we 
know that this is not so in the standard Yang-Mills case. 
Therefore, wc can already set = Q as this term is 
not allowed in the counterterm of the standard Yang- 
Mills action, which is a special case of the action we are 
studying[47]. Secondly, duo to the Landau gauge condi- 
tion (3.) and the antighost equation (4.) we find, 



ai = a2 



(A36) 



Next, the linearly broken constraints (5.) give the fol- 
lowing relations 

ai = -as = -0-9 = aio = nn = -63 = 64 , 

04 = as = — ae = aj , hi = 1)2 = = ~ . (A37) 

The Rij symmetry does not give any new information, 
while the integrated Ward identity relates the two previ- 
ous strings of parameters: 

ai = -as = -ag = aio = an = -bs = b^ 

= as = 04 = —as = ae . (A38) 

Taking all this information together, we obtain the fol- 
lowing counterterm 



renormalization of the fields and sources. If we try to ab- 
sorb the counterterm into the original action, we easily 
find. 



= 1-/^1 



zT = 1 + 



h(^+ai), (A39) 



and 



zlj' = zli- = z-^'i'z-^i- = 1 - /if , 
Zh = z^ , 
Zk = Zll\ 

z^ = zT. 



(A40) 



The results (A39) are already known from the renormal- 
ization of the original Yang-Mills action in the Landau 
gauge. Further, we also obtain 





= 4/^ = z-y^z-'/' = 1 




zT 


= z,^/\ 




z'J' 






Zm 






Zn 


= z-/^\ 




Zu 


= i + h'f = z;\ 




Zv 


= ^-h^i = z^''^z-/'\ 




Zj" 


= 1 + /^? = 




Zr 


= i-h'-^ = z-^i-z-^'\ 





(A41) 



This concludes the proof of the renormalizability of the 
action (AI) which is the physical limit of Sq^. 



S'^ = aoSyM + ai 



- gfabcR;%c''co^ + gfabcT^%c'Tp^^ . 

4. The renormalization factors 

As a final step, wc have to show that the countert- 
erm (B3) can be reabsorbed by means of a multiplicative 



Appendix B: Renormalization of the 
Gribov-Zwanziger action, option 2 

Wc can also write down an alternative proof of the 
renormalization of the Gribov-Zwanziger action. Look- 
ing at (20), we are tempted to treat gfauA'^ip'l^ and 
gfakbA^^Tp^J^ as the relevant composite operators. How- 
ever, we shall see that these operators mix with 9^<p"'^ 
and dfjjpl'^, respectively. A similar observation holds 
for gfakbA'i^ujl'' and gfakbA^^uj^". Such conclusion is ev- 
idently not unexpected, as we have learnt already from 
Appendix A that e.g. 

Dfipl" = d^iff + gfakbAf^ifi" is multipUcatively renor- 
malizable. 
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1. The stcirting action and the BRST 



2. The U{f) invariance is easily adapted 



Looking at (A5), we may try to also start with another 
possible BRST invariant action, 



(B5) 



^(2) 
^GZ 



Sym + 5gf + So + SP + S, 



ext ; 



(Bl) 



whereby wc have immediately added ^ext still given by 
(A9) and with 



j d'^x i-M. 

ahcjjfai j^bd 



akb A k b 



fc, h 



-gr^U'^W'^c^^'r + U'^-^gfakbA^oj. 



-T^^gUbcD'lfc^Tpf) 



(B2) 



In contrast with the previous section, we have now intro- 

duccd 5 doublets, Mp, {U';\ M';'), {V^\ Nf), 

X'"', N'^') and {T^\ Rf) with the following BRST 
transformations , 



^ 5M'>, 



^ su. 



+n:: 



I L yai 



v:: 



5C 



^ " 5R<f, 



^ ST: 



We have again that the i-valued fields and sources 
turn out to possess an additional quantum number. 
All the quantum number are still the same as in 
Tabic 1 and Tabic 11, whereby wc keep in mind 
that the quantum numbers of the primed sources 
are obviously the same as those of the unprimed 
ones. 

3. The Landau gauge condition does not change. 









= 0, 






sM'f 


= 0, 




— TV"' 




= 0, 




— TV'"* 




= 0, 










pai 

- ^11 . 




= . 



(B3) 



(2") 

In order to go back from Ss to Sg from Appendix A 
(sec equation (A6)), we just need to set U = U' , V = V', 
N = N' and M = M' . Eventually, it appears natural to 
give the primed sources the same physical value of their 
corresponding unprimed counterparts, see equation (A8). 



2. The Ward identities 

Just as in Appendix A, we enlist all the Ward identities 
obeyed by Sq2, which of course look very similar. 

1. The Slavnov- Taylor identity is now given by 

5(sg))=0, (B4) 



with 



GZ> 



GZ 



Sc" 



5K- 

.(2) 



5A- 



5c 



I —a '^^GZ I , ,0 



(2) 
GZ 



5^1 



R' 



.(2) 



ft 



I pr/ai "^GZ 



Shu - ^M^M 



(B6) 



4. The same goes for the antighost equation. 



Sc"- 



-a,— ^ = 



(B7) 



5. The linearly broken local constraints now become 



^^GZ , Q '^^GZ 



J-nf T"d» GZ 
-r gjdba-l-fj. ^^ftj 



9f 



abc Ab T//ci 



(B8) 



•^^GZ _|_^ <^^GZ 



6lj? 



SN°, 



-9f 



g <^^GZ 

,5E 



UJ 



(2) 

GZ _ 



= 5/ 



a6c^6 ^'ci 



(B9) 



We also find some extra linearly broken constraints 



(BIO) 
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6. The exact TZij symmetry can be adapted to 



5. The Unearly broken local constraints 



(Bll) 



with 



SN°, 



+V, 



uy 



u. 



5R7 



(B12) 



7. The integrated Ward identity is now linearly broken 
as follows 



5ujf 



rjlai " ^GZ 



(B13) 



5 

Wi 

6 



5 



5M?, 



gfabcT 



SK?, 



+5, 



_5_ 

'5¥ 



x4t = o 



and 



6. The exact 7^^ symmetry 



(B18) 





= 0, 




= 0, 






= 0, 




= . 


(B19) 



The counterterm 



We again translate all the identities into identities for 



the counterterm S, 



(2)c 
GZ 



1. The linearized Slavnov- Taylor identity: 



(B14) 



with S'^^ the nilpotent linearized Slavnov- Taylor 
operator, 



7? y(2)c 



. 



(B20) 



7. Finally, the integrated Ward identity becomes 



/ 



(2)c 
GZ 



(2)c 
GZ 



(2)c^ 



I rr/ai^^GZ 



. 

(B21) 



Now we can write down the most general counterterm 

^GZ*^ of (i = 4, which obeys the linearized Slavnov- Taylor 
identity, has ghost number zero, and vanishing Q f num- 
ber, 



AAA 

-|-Ar,f -— -h M'f—— + N'"' 



2. The U{f) invar iance 



^U^GZ 







ST" 



+ 



3. The Landau gauge condition 

= 0. 



^4t 



4. The antighost equation 



(B15) 



(B16) 



SK", 



(B17) 



Notice that the part in a and h parameters is exactly the 
same as in the previous Appendix A, see equation (ASS). 
We shall now impose all the constraints induced by the 
Ward identities. We keep in mind that the argument con- 
cerning the broken ghost Ward identity still holds. Also, 
the 4 constraints (B19) invoke the counterterm to be in- 
dependent of the sources [/', V , M' and N' . Ultimately, 
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we find 



Sq2^ = uoSym + ai 



check that equations (A39) and (A40) still hold, and also 
the Z-factors of zl'"^ , Z^"^ , Zt and Zr do 

not change. Only the renormalization of the sources U, 
V, M, N is different as they mix with respectively U', 
V\ M\ N'. Indeed, we find 'that 



(B22) 



We notice the close similarity between this counterterm 
and the one in expression (B3). 



The renormalization factors 



Mo 
Ml, 







-ai 
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The last step is to find all the renormalization factors. 
Due to the close similarity with the output of Appendix 
A, many Z factors will be the same. One can indeed 



which again proves the renormalizability of the Gribov- 
Zwanziger action. The consequences of this mixing shall 
be explained in sections II. B and III.B. 
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